arXiv:1507.06289vl [math.AP] 22Jul2015 


A FRACTIONAL FREE BOUNDARY PROBLEM RELATED TO A 

PLASMA PROBLEM 

MARK ALLEN 


Abstract. We study a fractional analogue of a plasma problem arising from 
physics. Specifically, for a fixed bounded domain we study solutions to the 
eigenfunction equation 

(-A) s u = A (u - 7)4. 

with u = 0 on 9ft. 


1. Introduction 

A mathematical model for the region inhabiated by plasma in a Tokamak ma¬ 
chine is given by the two dimensional equation 



in a bounded domain Q. u+ is the positive part of the function it. The region 
inhabited by the plasma is given by {it > 0}. Properties of solutions to this and 
similar problems was studied in m and [3]. The simplified model with solutions 
to 

(1.1) — A u = Au+ 

was studied in M- Regularity of the free boundary d{u > 0} was studied in 
m- The physical applications of the simplified model CCD are for two dimensions; 
however, one may study (ED in higher dimensions. The regularity of the free 
boundary for higher dimensions was studied in [10] . 

In this paper we study a fractional analogue of m- For a bounded domain f 1 
in M” we consider solutions to the equation 

(1.2) (—A) s it = X(u — 7)+ 

with it = 0 on dCl. The reason for the translation (it —7) is to utilize the definition of 
the fractional Laplacian with zero dirichlet data. In Section [2] we use the extensiion 
operator and subtract out the constant 7 which gives an equation of the sort 

“(-A ) s u = A u+ " 

The main aim of this paper is to study properties of the free boundary {it(-, 0) = 0}. 


M. Allen is supported by NSF grant DMS-1303632. 
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1.1. Outline. The outline of the paper is as follows: In Section [2] we establish 
certain properties of the fractional Laplacian that will be needed in our paper. We 
also discuss the notion of the extension operator that allows one to “localize” the 
fractional Laplacian. In Section[3]we prove existence of solutions to (11.21) . In Section 
[4] we prove interior regularity for solutions. In Section [5] we begin the study of the 
free boundary. We prove topological properties of the free boundary and show how 
they may differ from the original local plasma problem G3J. In Section [6] we use 
an Almgren’s type frequency function to classify so-called blow-up solutions. The 
classification of blow-ups allows us a classification of the free boundary points. We 
then give a regularity result for the “regular set” of the free boundary. In Section [7| 
we define the singular set and prove a Hausdorff dimensional bound for the singular 
set which shows that the singular set is “small”. 

1.2. Notation. The notation for this paper will be as follows. Throughout the 
paper 2s = 1 — a and — 1 < a < 1 and s will always refer to the order of the 
fractional Laplacian (—A) s . ( x,y ) £ R n+1 with x £ R" and y £ R. O will always 
be a smooth bounded domain. For a set U £ R” +1 , 

- L 2 (a, U) := {/ | f\y\ a/2 £ L 2 (U)}. 

- L 2 (a,dU) := {/ | /|y|“/ 2 £ L 2 (dU)} with respect to H n Hausdorff measure. 

- H^(a,TJ) := {/ | /, V/ £ L 2 (a, [/)}. 

- U' := {x £ R” | 0,0) £ U}. 

- We refer to R™ x {0} as the thin space. 

- B r := {x £ R n | 0| < 1} 

- f± denote the positive and negative parts of x respectively so that / = /+ — /_. 
We denote the free boundary as T = {u(-, 0) = 0}. 

2. Fractional Laplacian 

We define the fractional Laplacian through the spectral decomposition. For 
a bounded domain H let 0 < Ai < A 2 < ... and {</)} be the eigenvalues and 
corresponding orthonormalized eigenfunctions with dirichlet zero boundary data. 
For / £ L 2 (S2) we write 

(2.1) / = y^afc^fc- 

k=1 

Then the fractional Laplacian is given by 

( 2 . 2 ) (-Ar/oHE^^)- 

2=1 

The fractional Laplacian can also be given as a Dirichlet to Neumann boundary 
data map by the use of an extension operator. In the case when (—A) s is defined on 
all of 1" instead of on a bounded domain this equivalency was given in the paper 
[4j . For a bounded domain there is an analoguous extension operator 12]. We look 
at the solution to the following weighted elliptic problem in an extra dimension 

div(y a Vu) = 0 in fi x R + 
u(x, 0) = f(x) 

u(x, y) = 0 for (x, y) £ dfl x R + 
u(x, y) —» 0 as y —» 00 
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where a = 1 — 2s. The fractional Laplacian is a Dirichlet to Neumann boundary 
data map: 


(2.3) c a lim y a u y (x, y) = (~A) s u. 

y ->0 

where c a is a constant only depending on a and dimension n. If we shift the solution 
downward by subtractin the constant 7 our equation of study then becomes 

c a lim y~ a u v (x, y) = A u+ for every x £ fi. 
b-> 0 

Absorbing the constant into the right hand side we obtain 


(2.4) lim y a u y (x,y) = —A u+ for every x £ O, 

y-*- 0 

where A > 0 is a new constant. Many of the results in this paper will apply to 
solutions of the more general equation m- However, if u is a solution to Ol 
then after subtracting 7 , the extension function u will satisfy m, but will also 
satisfy the additional condition 


(2.5) u y < 0. 

Condition (12.51) will be used in Section [7] to give a Hausdorff dimensional bound on 
the singular set. That solutions of CuD satisfy (12.51) is easily seen from the following 
argument. If u is the extension for a solution to (11.21) . then y a u v is — a-harmonic 
(SI)- Furthermore, y a u y —> 0 as y —> 00 . Then y a u v has nonpositive boundary 
data on d(Ll x R + ), and so from the maximum principle for —a-harmonic functions 
( 0 ) we conclude that u y < 0 . 

We will utilize the following notion of trace for the weight y a (see EH)- 

Proposition 2.1. Let U C R n+1 an open Lipschitz domain. Then there exists two 
compact operators 

Ti : H\a,U) ^ L 2 (a,dU) 

T 2 : H\a,U) ^ L 2 (U') 


By utilizing rescaling and Proposition 12.II on B 1 we obtain the following 


Corollary 2.2. Let v £ H l (a,B+). Then there exists a constant C = C(n,a) such 
that 


f y a u 2 dTL n < Cr f y a \X7u\ 2 
J(dB r )+ JB+ 

[ u 2 dH n < Cr 1 - 11 [ y a \S7u \ 2 
JB' r J B+ 

We will also need the following Hopf type Lemma 


Lemma 2.3. Let v be a non-constant a-harmonic in B+ for some r > 0. Assume 
v achieves its minimimum at (xq , 0 ) £ B' r . Then 


lim 
v~* 0 


’{x 0 ,y) 

n.l— a 


> 0. 


Proof. By subtracting a constant we may assume v(xo, 0) = 0 and therefore v > 0. 
Let w (not identically zero) be an a-harmonic function satisfying w(x, 0) = 0 and 
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w < v on dB r . By the Boundary Harnack Principle for a-harmonic functions stated 
in [5; we have for p <dist(:ro, dB r )/2 


w 


Then 


and so 


sup —:— < C(r) inf ——. 

B p (x 0 , 0 ) y “ B p (x o ,0) y a 


i»f rfe»)>o, 

b p (x 0 , o) y l ~ a 

lim A™! > lim > o. 

y—^ 0 y 1 a y—>-0 y 1 a 


□ 


Remark 2.4. The Boundary Harnack Priniple is a powerful tool, but is not necessary 
to prove Lemma 12.31 One can for instance use the power series representation for 
odd a-harmonic functions shown in [ 1 : to achieve the same result for w in the above 
proof and hence for v. 


3. Existence 


In this section we prove existence of solutions to m- To obtain an eigenfunction 


we consider minimizing the 

fractional energy 


(3.1) 

D{u) := f u(- 

-A ) s u 


Jn 


subject to the constraint 



(3.2) 

G(u) := f (u- 
Jn 

7)+ = 


where c, 7 > 0 are two fixed constants. Using the extension mentioned in Section 
[21 this is equivalent to minimizing 

JJ y a l Vu i 2 dx d l b 

f2xR+ 

subject to the constraint (13.21) . 

Lemma 3.1. There exists a minimizer of G3D subject to the constraint (13.21) . 


Proof. Minimizing is equivalent to minimizing 

JJ y a \\7w\ 2 dx dy 


QxM+ 


with w(x,Q) = u(x) and w(x,y) —> 0 as y —>• 00 . From Corollarv l2.2l it follows that 
L 2 (U) is compactly contained in H 1 (a , fl x R + ). The existence of a minimizer then 
immediately follows. □ 


We note that the extensions is not necessary to prove a compactness theorem. 
Using only the spectral decomposition there is an elementary proof using power 
series that if a sequence Uk, is bounded in then there exists uq and a 

subsequence such that Uk — 1 uo in H s (£l) and Uk —> uo in L 2 (f 2 ). 

Lemma 3.2. For a domain U and fixed constant c there exists a solution to (L2} 
for some A > 0. 
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Proof. Since the functionals D : H s —> R and G : H s / 2 —> R have the Frechet 
derivatives 

D\u) = (-A ) s / 2 u G'(u) := 1/2 [ (u- 7 )+ 

Jn 

then there exists A > 0 such that for a minimizer of D(u) subject to the constraint 
G(u) = c satisfies 

D'(u) = XG'(u). 

See [15] for a discussion of Lagrange Multipliers with Frechet derivatives. □ 

The above Lemma shows existence for some A > 0. In m existence is shown for 
fixed A belonging to a correct range. It would be of interest to show existence for 
fixed A for this fractional problem as well. Since this paper focuses on properties 
of the free boundary, we chose to only give a quick proof for some A to simply 
demonstrate that our class of solutions we study is nonempty. 

4. Interior Regularity of Solutions 

In this section we obtain the interior regularity of solutions which will enable us 
to obtain regularity of the free boundary where the gradient does not vanish. Since 
we are dealing with an eigenvalue equation, we use a bootstrap technique. Obtain¬ 
ing regularity for u passes the same regularity to A u + up to Lipschitz regularity. 
Regularity for Au + = “(— A) s u" allows us to obtain higher regularity for u. To use 
the bootstrap technique we utilize the following Proposition from [Ill- 

Proposition 4.1. Let w = (— A) s u. Assume w £ C 0 , “(R”) and u £ L°° for 
a £ (0,1) and s > 0. Then 

If a + 2s < l then u £ C°’“ + 2 s (R"), 

If a + 2s > 1 then u £ C 1 ’“ + 2 s- 1 (R ra ). 

Using Proposition 14.11 we are able to obtain the following interior regularity 
result. 

Theorem 4.2. Let u be a solution to (|2.4D . Then for any I\ (e fl 

If s < 1/2, then u £ C 1,a (K) for every a < 2s, 

If s > 1/2, then u £ C 2,a (K ) for every a < 2s — 1 

Proof. Throughout the proof K will be any compact subset of Q. Each time we 
obtain higher regularity in the bootstrap technique we will pass to a smaller subdo¬ 
main. By abuse of notation we will continue to call the smaller set K. The Holder 
norms will of course depend on I\. To utilize Proposition l4.ll which applies to the 
fractional Laplacian whose domain is functions defined on all of R”, we use the 
Reisz potential as well as the a-harmonic extension of u in R ra x R + . 

u 0 ( Z ) := (-A)-/A u+ = \ z X ™J%l 2a dy- 

If we consider the variable 2 £ 1" x R + , then uq is a-harmonic in l n x K + . 
Furthermore 

lim y a d y {u — uo)(x, y) = 0 for any x £ K 
2/->0 

By even reflection in the y variable u — uq is a-harmonic in the cylinder K x R. 
a-harmonic functions have a power series representation and hence are smooth [Jj. 
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Therefore, since u — uq is smooth on any K <g K , any regularity obtained for wo 
on K is passed to w on K. Initially we only know that u + G L 2 (Lt). From the 
theory of Reisz potentials, wo G L q (I\ ) for any 1/2 — 1/q < 2 s/n, son £ L q (K), 
and so w + G L q (K ), and we rename K as K and continue the process. Then using 
again imbeddings of Reisz potentials we obtain that after finitely many iterations 
(depending on s), « £ L°°(K) and hence w + as well. Then from [TT] we obtain that 

w 0 £ C°’“(R n ) if a < 2 s < 1 

w 0 G C , 1 ’“(R") if a < 2s - 1 and 2s > 1 

Again, we obtain the same regularity for w inside K. Having obtained initial Holder 
regularity for w + , we obtain higher Holder regularity for u by using Proposition l4.ll 
After finitely many iterations (again depending on s) we obtain the conclusion of 
the theorem. Lipschitz regularity is the most we can know for w + , and this is when 
the iteration process stops. □ 

5. Topology of the Free Boundary 
For the original local plasma problem studied in HD], the free boundary is exactly 
d{u > 0 } = <9{w < 0 } = {u = 0 }. 

That these two boundaries are the same follows from the maximum and minimum 
principle since Aw = 0 in {w < 0} and Aw < 0 in {w > 0}. The situation is different 
in the nonlocal/fractional case: we cannot apply the same minimum principle to 
solutions of 

lim y a u v {x,y) < 0 
v-y o 

since it is possible to have a local minimum and still satisfy the above inequality. 
It may be possible to construct solutions to (ED that satisfy 

d{w( • ,0) < 0} C d{u{ • ,0) > 0}. 

We are mostly interested in the portion of the free boundary d{u( ■ , 0) < 0}. This 
next proposition gives an inclusion when we assume additionally ESJ). 

Proposition 5.1. Let u be a solution to (12.411 and assume also (12.51) . then d{u < 
0} C d{u >0}. In particular, if u is a solution to (11.21) . then d{u < 7 } C d{u > 7 }. 

Proof. Let V be an open subset of LI such that V C {w( ■ , 0) < 0}. Suppose there 
exists xq G V such that (xo,0) G 9{w(-,0) < 0}, and so u(xq,0) = 0 and w is not 
constant. Now y a u v (x,y ) = 0 for x G V. If we evenly reflect in the y -variable u 
is a-harmonic on the domain U x (—e, e) for e > 0. From E3 it follows that u 
achieves an interior maximum at (cco, 0 ), and since w is not constant we obtain a 
contradiction to the maximum principle for a-harmonic functions [ 8 |. □ 

This next Proposition shows that if w is a solution to UE2D. then w is strictly 
subharmonic - in the classical sense for the Laplacian (not fractional Laplacian) on 
the thin space - across the boundary d{u > 7 } which is not true when s = 1 . 
This illustrates why one cannot hope for a strong minimum principle in {w > 7 }. 

Proposition 5.2. Let u be a solution to E3 with 1/2 < s < 1. Then there exists 
an open set U (in the topology of LI) containing |w(-,0) = 7 } such that 

A.u(x, 0) > 0 for every x G U, 

where A is the classical n-dimensional Laplacian on the thin space R”. 
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Proof. Since s > 1/2, from Theorem 14. 2l we know u is C 2 on the thin space R", so 
Au(-, 0) exists in the classical sense. Also, from the extension we have that 

v(x,y) = c a y a u y (x,y) > 0, 

where c a < 0 is as in (12.311 . and v is — a-harmonic. Let u(X q, 0) = 7 . Then 

lim y a u y (x 0 ,y) = 0 . 

V ~*0 

We then conclude that (xo,0) is a minimum for the —a-harmonic function v. and 
so from Lemma 12.31 

0 < lim y~ a v v (xo,y) 

y ->0 

= -(-A) 1 -'«(a:o) 

C a 

= — (—A) 1 _s (—A) s u(xo) 

C a 

= —(-A )u(x 0 ). 

Ca 

Since c a < 0, we conclude that Au(x, 0) > 0 whenever u(x, 0) = 0. From the 
C 2 continuity of u it follows that u is strictly subharmonic in a neighborhood of 
{w = 7 }. □ 

In the next theorem, we show that if Q has symmetry in a coordinate direction, 
then the same symmetry will be inherited in the level sets of solutions of ca that 
arise as minimizers of EED subject to the constraint (13.21) . 

Theorem 5.3. Let u be a solution to OD that arises as a minimizer of EED 
subject to the constraint m - Assume Q is symmetric across a hyperplane P 
orthogonal to the direction v. Then the level sets of u are also symmetric with 
respect to P. 

Proof. We pick our coordinates so that P = (0, X 2 , ■ ■ ■, x n ). We use the a-harmonic 
extension. Let u minimize the energy functional 

/ v a |Vu| 2 

JQxR 

subject to the constraint (13.21) . If we Steiner symmetrize in the X\ direction 
(which is orthogonal to y) to obtain u* , and if 0 < Tj < T 2 < 00 , then 

(5.1) II y a \Vu* | 2 < II y a \Vu \ 2 

fix[Ti,T 2 ] nx[T!,T 2 ] 

and equality is achieved if and only if u is already Steiner symmetric in the x\ di¬ 
rection. The case of equality holds from the result in [9] since a-harmonic functions 
are real analytic away from y = 0. That is why initially we restrict ourselves to 
Xj > 0. Now we let Tj —► 0 and Xj —>• 00 to obtain 

II y a \S7u*\ 2 < II y a \S7u\ 2 

OxR+ JIxl+ 

If u is not already Steiner symmetric in the Xi direction when y = 0, then by 
continuity u will not be Steiner symmetric at some time X > 0, and so by (ED 
the energy of u* will be less on an interval x |Tj, X 2 ] and hence also on x R + . 
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Notice that the constraint (13.21) is preserved for u*. Then if u is a minimizer of 
(EH), it must be symmetric in the xi direction. □ 


This next Corollary gives a sufficient condition on the shape of the domain f l 
under which d{u > 7} = d{u < 7} for solutions of ( 11 . 21 ) that arise as minimizers. 
A good question would be what conditions are necessary on fl in order to assure 
the same condition on the free boundary. 


Corollary 5.4. Let u be as in Theorem 1 5. 31 and assume f 1 is symmetric with 
respect to 27 for 1 < i < n. Then 

d{u( • , 0) < 7} = d{u( • , 0) > 7}. 


Proof. From Theorem 15.31 u is symmetric in each 37 variable. Thus u achieves a 
maximum at the origin. If u(0) < 7, then u is a-harmonic everywhere and this is 
a contradiction to the maximum principle. Then u(0) > 7. Suppose there exists 
z = (zi,Z2, ■ ■ ■ , z„, 0) € d{u(-, 0) > 7} and 2 ^ d{u( 0) < 7}. By the symmetry of 
u we can assume 27 > 0 for each i. Also from the symmetry of u we have 



< 0 . 


For each x in the first quadrant and in particular for the point z. If z d{u < 7}, 
then by continuity there exists an open region 


U := {x \ x £ B e (z ) D {27 > Zi} } 

and zt(x, 0 ) = 7 for every x G U. u — 7 is then an a-harmonic function with 
both an even and odd extension (in the y-variable) across the thin space. Then 
from the power series representation of a-harmonic functions (see m) it follows 
that u is constant. This is a contradiction since u(x,y) -A 0 as y —>• 00 , and 
u(0, 0) > 7 > 0. □ 


6 . Regularity of the Free Boundary 


In this section we look at the regularity of the free boundary. Our results apply 
to solutions of (12.41) . We now state a Lemma that will allow us to utilize Almgren’s 
frequency function. For solutions of (12.41) Almgren’s frequency function will not be 
monotone. However, we will prove that the limit at the origin exists and use this 
result to put a bound on the dimension of the singular set of free boundary points. 
We define 

D(r) := [ y a |Vw | 2 H(r) := f y a u 2 N(r) = r 

Jb+ J(dB r )+ H(r) 

Lemma 6.1. Let u be a solution to (12.41) . Then lim r _ >0 + N(r) exists and 

0 < lim N(r) <00 

r —>0 


Proof. For solutions of (12.41) . N(r) will not necessarily be monotone. We therefore 
begin by considering the modified function 


N(r) 


d( t) - a Ib' u 2 

r H(r) " 


Both N(r), N(r) are absolutely continuous and hence differentiable for almost every 
r. We note that N(r) > 0. By the same computations as in [J and using that u 
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is a solution to (12.411 (along with the accompanying C 1,a regularity) we obtain the 
following Rellich-type identity 

(6.1) D'(r) _ n ~ 1 + a D(r) + 2 [ y a u 2 v + - f (x,Vu 2 + ). 

r J(dB r )+ r J B , r 

We also have by routine computations 


,, . n + a , , f 
H (r) = ——— H{r) +2 y 

r J(dB r )+ 


( 6 . 2 ) 


d_ 

dr 


’OB' 


[ (a:,V4) = r f u\ - n [ 
J B' J dBL J E 


Jl +- 


Combining EH) and m we have 

d , ~ . . 1 n — 1 + a 

— log TV r = - + 

dr r 


D(r) 


r D{r) — Xf B , 


b: “+ 


+ o I(dB r )+ y a ' u l _ 2 f(dB r )+ y‘ 


n + a 
r 

1 uu l , 


‘ D{r) + A f B , u\ f (dBr) + V a u 2 


nX 


Ib' u + 


r D(r) + \f B , u 2 + 


Using integration by parts and that u is a solution to (12.41) we have 


D(r) - A [ u 2 + = f 
Jbl ju 


Then 


d , . 1 — a 

— log N(r) = - 

dr r 


+ 2 


1 - 

'/I 


(9B r )+ 

D(r) 


y uu v . 


D ( r ) + Xf B , u\ 


(OB r ) 


+ y a ul 


I, 


(OB r ) 


+ y uu u 


Ju 


f(ffR.)+ y au2 


l(OB r )+y UUv J(OB r y t 

The first term is clearly nonnegative, and he second term is nonnegative by the 
Cauchy-Schwarz inequality. Since N(r) is monotone, the limit as r —> 0 exists, and 


0 < lim N(r) < oo. 

r^r 0 + 


(6.3) 

Now from Corollary 12.21 we have 

(1 - CAr 1_a )iV(r) < N(r) < N{r), 

and so (IQ) holds with N(r ) replaced by N(r). To show that the limit is greater 
than zero we note the following rescaling property N(r,u) = N(l,u r ) where 

u{rx , ry) 


(6.4) 


<i a oi 2 ^ 


1 / 2 - 


(j^f(0B r) + y“ u 

From the rescaling and since N (r) is bounded we have 

I H 1 (a,B+) — 


Ur\\L^(a,(0B 1 )+) = 1 ||W- 1 | pi f n R+^ < M 
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where M is a constant depending on u. Then from Proposition 12.II we have that 
there exists a sequence rk —> 0 and a function uq such that 


u rk u o in H l (a,B+) and u Tk —> uq in L 2 (a, ( dB r ) + ). 


We note also that 


Now 



2/°|Vu 0 | 2 


N( 0+). 


[ y a u o 

J(dB r )+ 


= 1 , 


so uq is not identically zero, and hence we conclude that N(0+) > 0. 


□ 


Corollary 6 . 2 . Let u be a solution to (12.41) . Let u r be defined as in (16.41) . Then 
for every sequence rk —> 0 , there exists a subsequence and a function uq so that 

u Tk — 1 uq in H 1 (a, B+). 

Furthermore, N(u, 0+) = N(uo,r) for every 0 < r < 1 and if we evenly reflect uq 
in y, then u o is a-harmonic and homogeneous of degree N(u, 0+). 

Proof. The weak convergence of a subsequence was established in the proof of 
Lemma l 6 Jl Now for t > 0 


N(t,uo) = lim N(t,u r ) = lim N(tr, u) = N(u, 0 +), 

r—>-0 r—>-0 

so N(r,uo) = N( 0+). From the rescaling we have the property 

lim y a d y u r (x,y) = -rX(u r )+. 
v-> o 

If we let ip £ Cq(Bi), then 


hm y a (Vip,X7u rk ) =r k \ (u rk )+ip 0 . 

r ^°J B + Jb[ 

The right hand side goes to zero by utilizing Proposition 12.II Then if we extend u 
evenly in the y variable, u is a-harmonic in B\. □ 


The blow-up solutions will satisfy ( 12 . 41 ) with A = 0 . We can therefore evenly 
reflect across the thin space {y = 0}, and each blow-up solution will be even and 
homogeneous. From [6] it follows that the blow-up solutions will be homogeneous 
of order k £ N. When we have the additional assumption u y < 0 - as is the case 
when u is a solution to m - then we have a better classification. 

Corollary 6.3. Let u be a solution to ( 12 . 41 ) and let Xq £ T. Assume also u y < 0 . 
Let u r -A v be a blow-up of u at xq. If V x u(xq) ^ 0 , then v is a linear function. If 
S7 X u{x o) = 0 we have the following alternative, either (/): v is either homogeneous 
of degree 2 and of the form 

(6.5) v{x, y) = p(x) — cy 2 

where p{x) is homogeneous of degree 2 in the x-variable and c > 0, or (II): v y = 0 
and A«(•, 0 ) = 0 and u is homogeneous of degree k £ N with k > 2 . 
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Proof. When V x u(xo) 7 ^ 0 this is immediately clear since v = V x u(xo). Now 
since y a u y (x,y) < 0 for all (x,y), this inequality is preserved in the limit for v. 
Then y a v y is —a-harmonic, nonpositive for y > 0, and has zero dirichlet data when 
y = 0 (since v is even). From the Boundary Harnack Principle 0 it follows that 
y a v y is comparable to y 1+a , or identically zero. If v y = 0 then the classification is 
immediate. If v y is not equivalently zero, then since v is homogeneous y a v y is also 
homogeneous, and thus it is homogeneous of degree 1 + a. Then v is homogeneous 
of degree 2. Since v y < 0 for y > 0 it follows that v must be of the form (16.51) . (See 
also p] for a classification of homogeneous a-harmonic functions). □ 

We now define the regular set of the free boundary. Let u be a solution to (12.41) . 
R u := {x | u(x, 0) = 0 and X7 x u(x,0) 7 ^ 0}. 

From Corollary [672] R u is equivalently the set of free boundary points where 
N(u, x, 0+) = 1 or equivalently the set of free boundary points that has a blow-up 
(and hence every blow-up) that is a linear function independent of the variabley y. 

From the implicit function theorem and Theorem 14.21 we have the following 
regularity result for R(u). 

Theorem 6.4. Let u be a solution to flU) . If u(x 0 , 0) = 0 and (x-q, 0) £ R u , then 
in a neighborhood V of x 0 , {u(x , 0) = 0} is a C 1,a (C 2,a ) graph for s < 1/2 (> 1 / 2 ) 
and a as in Theorem ]^. 2\ 

It is common in free boundary problems for the free boundary to be more regular 
than the solution. For instance around regular points for the original local plasma 
problem m the free boundary is real analytic [TO]. A question of interest would 
be to prove higher regularity of the free boundary for solutions of (12.41) . 

7. The Singular Set 

In this section we consider solutions to (El with the additional assumption 
u y < 0. This is the case for instance when it is a solution to El- The main result 
of this section is to give a Hausdorff dimensional bound for the singular set defined 
as 

:= {(*,0) £ 3{u(-,0) < 0} | S7 x u(x, 0) = 0}. 

Notice that we do not consider the set (x, 0) £ (9{it(-,0) > 0} \ 9{u(-,0) < 0}). 
The main result of this section is 

Theorem 7.1. Let u be a solution to m in Q. Then the Hausdorff dimension 

H T (S U ) = 0 

for t > n — 2. 

To prove Theorem l7.1l we follow Federer’s method of dimension reduction utilized 
in [7j. The Theorem is an immediate consequence of the following two Lemmas. 

Lemma 7.2. Assume that 'H T (S V ) = 0 for all functions v which arise as blow-ups 
of solutions to (12.41) . Then TL T {S U ) = 0 for all solutions u to (12.41) in ft. 

Proof. We first show the following Property (P): for every x £ S u there exists 
d x > 0 such that for all 5 < d x any subset D of S u D Bg{x) can be covered by a 
finite number of balls B ri {xf) with Xi £ D such that 
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We show (P) using a compactness argument. Suppose by way of contradiction the 
conclusion is not true for a sequence 5 k —>• 0. By picking a subsequence if necessary 
we perform a blow-up at x, us k —> ug. By assumption there exists finitely many 
B ri / 4 {xi) such that 

S Uo C UP n/4 and < -. 

i 

Since u< 5 fc — » ug in C 1 ’^ it follows that there exists kg such that if k > kg then 


s USk n Pi c u B ri/i . 


Then rescaling backward u satisfies the conclusion for all large k in B$ k and we 
reach a contradiction. 

We now denote D k := {y £ S u \ d x > 1 /k}. Fix x a £ D k . By Property (P) we 
can cover D k fl B ro (xg) where r 0 = 1/k with a finite number of balls B ri (xi) with 
Xi £ D k and 


E 



Now we repeat the same argument for each B ri {xi ) and cover it with balls B ri . (x^) 
to obtain 


E 



3 

Repeating this argument m times we obtain 'H s {Dk H P ro (xo)) = 0. 
k —> oo to conclude the Lemma. 


We then let 
□ 


Lemma 7.3. Let ug be a blow-up of u k {x o) where u is a solution to (12.411 and 
x Q £ S u . Then 

TL T (S Uo ) = 0 for t > n — 2. 

Proof. We begin by considering the problem in dimension n = 1, and suppose there 
exists Xo £ S u . From Corollary I6.3l un = ax 2 — cy 2 . Since Xo £ <9{u(x,0) < 0}, it 
follows that 0 £ 9{rto(a ; I 0) < 0}. Since c < 0 we have a > 0 and we immediately 
obtain a contradiction. The result now follows by using the standard argument 
of Federer for homogeneous solutions. In dimension n = 2 if 0 ^ Xo £ S Uo then 
performing a blow-up at Xg and using the homogeneity of ug we obtain a blow-up 
limit u oo in n = 1 with 0 £ S Uoo which is a contradiction. The argument for higher 
dimensions then follows using the ideas of Lemma 17.21 as in [7] . □ 


We remark here that using the same arguments it is easy to show W (T) = 0 for 
r > n — 1, so that the free boundary is of Hausdorff dimension n — 1. 
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